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Abstract 

Multistep Bose-Einstein condensation of an ideal Bose gas in anisotropic harmonic 
atom traps is studied. In the presence of strong anisotropy realized by the different trap 
frequency in each direction, finite size effect dictates a series of dimensional crossovers 
into lower-dimensional excitations. Two-step condensation and the dynamical reduc- 
tion of the effective dimension can appear in three separate steps. When the multistep 
^ behavior occurs, the occupation number of atoms excited in each dimension is shown 

to behave similarly as a function of the temperature. Multistep behaviors can be easily 
controlled by changing the degree of anisotropy. 

1 Introduction 

Although Bose-Einstein condensation (BEC) has been an active topic of research in con- 
densed matter physics for decades [p], 0, |3|, |J, recent revival of interests in this field is mainly 
credited to the achievement in atom molecular optics; most recent laser cooling and evap- 
orative cooling techniques in magnetic and optical atom traps enable us to realize BEC of 
a weakly interacting gas in a controlled way [||, |B|, [/J. The weakly interacting nature of a 
gas allows us to handle the phenomena and make theoretical predictions with high accuracy. 
The rapid progress in this field also stimulates other areas in physics such as high energy 
physics and astrophysics ||. 
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For atom trap experiments, since atoms are trapped in a finite geometry, finite size as well 
as finite number effects play a significant role in the condensation process. The conventional 
phase transition picture defined in the thermodynamic limit has to be reexamined or modified 
H pi], [II], [12], [13], [14], [15], [16| . One significant change due to finiteness is the existence of BEC 
in low- dimensional systems. Recently, quasi low- dimensional systems prepared by optical 
or magnetic trapping devices are actively discussed [O, |I§|. The study of such systems 



provides an ideal test for the theory of finite size, low-dimensional systems in a controlled 
environment. 

The critical behavior of a finite size system fl9| , 



infrared dimension (EIRD) of the system [22] , |23 



21] is characterized by the effective 
near the critical point when the contri- 
bution of the lowest mode of a system dominates, its symmetry properties can be shown to 
be equivalent to a lower-dimensional one. The system, in such a case, is said to possess an 
EIRD. Varying the relative size (or shape) of the potential changes the infrared behavior 
and hence the EIRD of the system. 

Dimensionless parameters 7/j = (3TiuJi {(3 = 1/ksT) for a harmonic oscillator potential 
with natural frequencies u>i (i = 1,2,3) characterize the degree of anisotropy and finite size 
effects. With respect to r/j, we can classify the dynamical behavior of the system into the 
following four cases dependent on the degree of anisotropy: 



Case 1 
Case 2 
Case 3 
Case 4 



r ]1 ,r ]2 ,r]z>l - EIRD = 0, 
771,772 > 1 >Vs -> EIRD = 1, 
r h >l>r ]2 ,r ]3 ^EIKD = 2, 
l>Vi,V2,m - EIRD = 3. 



As the temperature is lowered, dynamical dimensional reduction of the system characterized 
by the decrease of effective dimension can be observed. Particularly, in the presence of 
maximal anisotropy 1 >> 771 >> ^2 » V3i EIRD decreases one by one from three to zero 
at the temperature k B T = fko\, hu)2, and hu> 3 . 

While the picture described above is of generic nature of any quantum system, for par- 
ticles that obey Bose statistics, the dimension of the dynamics of excited particles can also 
change in the form of condensation. BEC can also occur in separate steps in the presence of 
anisotropy reducing the effective dimension attributed to the excited modes of the system 
by one, two, or three at a time Jjl|. In this sense, the dimensional crossover associated 
with both the multistep condensation and the reduction of effective dimension due to frozen 
degrees of freedom can manifest themselves in a similar way in a finite system in spite of 
their different origin. This is the main subject discussed in this paper. 

For the system with finite size and number of atoms, the reduced chemical potential <fi = 
(3(E — /i) does not strictly vanish at the critical temperature. Only in the thermodynamic 
limit, <p vanishes at the critical temperature and the specific heat develops the discontinuity in 



the derivative at the critical point |20|. In an isotropic system, the thermodynamic limit can 
be uniquely defined as discussed in many textbooks ||. However, if we allow anisotropy in the 
system, there are three different ways of taking a thermodynamic limit: (Three-dimensional 
limit) u>i, uj 2 , C1J3 — > while keeping NU1U2OO3 fixed. In this case, the system shows the critical 
behavior of three dimension and the corresponding three dimensional critical temperature 



T 3 o can be denned. (Two-dimensional limit) u 2 u 3 - ► while No) 2 o> 3 fixed. The system 
shows the critical behavior of two dimension and the corresponding critical temperature T 2 d 
is the one for the two-dimensional system. (One-dimensional limit) If we simply take u 3 — > 
while Nu>3 fixed, the critical temperature vanishes. However, if we tune uj 3 a little slower 
such that co>3 ~ log(2N)/N — ► 0, the one-dimensional critical temperature T\d can still be 
defined in this modified sense. We will discuss this issue again in Sec. 3.1. 

In the presence of strong anisotropy, the whole particle spectrum naturally splits into 
zero, one, two, and three-dimensional excitations. The ground state is viewed as a zero- 
dimensional excitation. Let us denote the number of modes excited in the corresponding 
directions as A*o, N\, N 2 , and N3, respectively. An n-dimensional condensation temperature 
T n o (n = 1, 2, 3) can be defined as the temperature at which all the n-dimensionally excited 
modes are saturated: 

3-dimensional; A" = N 3 (T 3D ), (1) 

2-dimensional; N = N 3 (T 2D ) + N 2 (T 2D ) , (2) 

1-dimensional; AT = N 3 (T 1D ) + N 2 (T 1D ) + N^Tm). (3) 

One can see that the condensation temperatures defined above are equivalent to the critical 
temperatures if the appropriate n-dimensional thermodynamic limit is taken. This splitting 
of the excitation spectrum gives the basis of the rest of our analysis. Similar splitting was 
proposed in [H| for liquid helium. In a liquid, however, an occupation number of particles 



excited in a particular direction is extremely difficult to observe. Furthermore, the validity 
of this splitting is by no means obvious for the strongly interacting system such as the liquid 
helium. On the other hand, such a quantity is directly observable in atom trap experiments 
and therefore it deserves careful study. Moreover, as shown in Section 3, occupation numbers 
Ni, N 2 , and A3 behave as if they were independent quantities and show the similar behavior 
when the multistep crossover occurs. This result indicates the independent nature of each 
A^j for strongly anisotropic systems. 

For a realistic system where N, u>i, 0)2,0)3, are all finite and fixed, physically observable 
temperature of interest is the crossover temperature at which the deviation from the bulk 
critical behavior sets in. The crossover temperature is achieved when the correlation length 
reaches the size of the system since further ordering in this direction will be suppressed at this 
point. In the strongly anisotropic system in which 0)1,012 >> o> 3 holds, Tir>, T 2 d « T 3 o gives 
the necessary (but not sufficient) condition for the multistep condensation: the condensation 
into two, one-dimensional modes and into the ground state can occur in separate steps. The 
multistep condensation was discussed for a nonrelativistic ideal gas in a cavity [I0[j, in a 



harmonic trap [24], and a relativistic ideal gas in a cavity [f25|j . 

In Section 2, we study the excitation spectrum of anisotropic harmonic oscillators. We 
focus our attention to three different cases in which the equipotential surface has a prolate, 
oblate, and maximally anisotropic ellipsoidal shape. In Section 3, we show that each case 
shows qualitatively different condensation behavior. After introducing the condensation 
temperatures defined in the bulk limit, we focus on the multistep crossover behavior of 
excited modes between different effective dimensions through BEC or dynamical reduction 



of EIRD. In particular, in a maximally anisotropic potential, the dimensional reduction can 
occur in three-steps. In such a case, we show that each dimensional component behaves in 
a similar manner as a function of the reduced temperature defined near the corresponding 
crossover temperature. The possible effect of interactions is also discussed. 

This work deals with a nonrelativistic ideal Bose gas in anisotropic magnetic traps, and 



a companion paper deals with a relativistic ideal Bose gas in rectangular cavities ||25|| . Our 
calculations are focused on the occupation number and condensation temperature for each 
dimension. The effect of an interaction in condensation process has been discussed in many 
literatures. In principle, it can affect the dynamics of condensation considerably. For a 
weakly interacting gas, however, the averaged quantities such as condensation fractions and 
critical temperatures are relatively insensitive to the presence of interactions and the cor- 
rections to bulk ideal-gas value are well-explained by the finite number correction |26|, [2Tf . 
On the other hand, interaction effects are known to affect higher moments such as specific 
heat significantly and considered to be essential to explain the observed specific heat data. 
Throughout the rest of the paper, we use units such that ks = % = 1 for brevity. The results 
in ordinary units can be easily reproduced by replacing u — * Tioj and T — * ksT . 

2 Anisotropic Harmonic Oscillator and Excitation Spec- 
trum 

For an anisotropic harmonic oscillator with oscillator frequencies u>i (i = 1, 2, 3), the Hamil- 
tonian has the form 

H= 1 -J2(P*+^)- (4) 

/ i=i 

In this paper, we study cases where the frequencies tu, are rationally related.^] Whence there 
exist integers ki (i = 1,2,3) such that Uiki — u> (i — 1,2,3) where u = ^(k^ks) 1 ^ 3 and 
U>lU2U>3 = O 3 . 

The energy level of an anisotropic harmonic oscillator is given by 

3 

E n = Y J ^i{n i + 1/2). (5) 

We can also define the energy level modulo fc« as rii = kiVi + \ where z/j = [rii/ki], and [ 
denotes the integer part of the number inside the bracket. Then Eq. ([5]) can be written as 

3 

E n = uM + Y,"i>* + Eo, (6) 

where M = v\ + V2 + z/3. The first term in Eq. (|6|) corresponds to the isotropic harmonic 
oscillator Hamiltonian (see also Appendix A). The ground state energy E has the familiar 

1 This assumption of rationality is rather for the technical convenience and will not affect the physical 
results of this paper. Other methods can be found, for example, in urn, p3, and p8[. 



form 

^1+^2+^3 

E = • (7) 

2.1 Excitation spectrum 

2.1.1 Prolate shape potential 

First we discuss the case of anisotropy corresponding to i0\ = uo 2 > oj 3 (we simply choose 
here k\ = k 2 = 1 < k 3 ). In such a case, the equipotential surface has a prolate shape. For a 
strong anisotropy k 3 >> 1, two-step condensation can occur. In such a case, u = u)\ = oj 2 
and the energy eigenvalue is 

E n = uM + W3A3 + E Q . (8) 

For sufficiently large values of ^3, the whole energy spectrum can be split into the energy 
level of the ground state ( E n — 0), one-dimensionally excited states ( E n = n 3 u 3 ; n 3 = 
1, 2, • • •), two-dimensionally excited states ( E n = n\W\ + n 3 u 3 ; ri\ = 1, 2, • • • , n 3 = 0, 1, • • • 
and n 2 uj 2 + ^3^3, n 2 = 1,2, •• -,713 = 0, l,---), and three-dimensionally excited states ( 
E n = n 1 uj 1 + n 2 LU 2 + n 3 uj 3 ; ni,n 2 = 1, 2, • • ■ ,n 3 = 0, 1, ■ ■ • ), where E n = E n -E = luM+lu 3 X 3 
is the energy measured from the ground state. The number of particles excited in these 
dimensions are given respectively by 

N = ^-, (9) 

1 — z 
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n 2 = y - 

ni=l,n 3 =0 c A 

fc 3-! °° 2Mz 

= 1^ 1^ p M m +\ 3 rj 3 _ ' ( ll ) 

A 3 =0M=l e * 

n 3 = y — —A 

ni=n2=l,ri3=0 

^ « (M-l)M 

A 3 =0M=2 ^ c xi 

where z = eP^~ E{ ^ is the (reduced) fugacity The factor 2 in Eq. ([TT|) accounts for the 
symmetry between the first and second axis. These expressions can be further simplified in 
the following manner. 



-n 3 r?3 00 | — /J73 



For one dimension, 

oc 

Nl = y — = y *— (13) 

n 3 =l (=1 

= ^ e "" /2) + ..., (14) 

V3 



where g p (z) = Y^=i f^ is the Bose-Einstein function. We used Eq. ( |B4|) to obtain the second 
line from the first line. 



For two-dimensional excitations, making use of Eq. (|B"2j), Eq. (11) can be written as 



k$ — 1 oo oo 

N 2 = 2 J2 E E Mz l e- l ^+ x ^) 

A 3 =0M=1 1=1 
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To obtain the third line from the second line, we used Eq. ( |B5| ) . 
For three dimensional excitations, Eq. ( |T2"| ) gives 

JV 3 = V V V ( M - 1 ) M J c ^(Mr ?1 +A 3 r ?3 ) 
A 3 =0 M=2 «=1 ^ 

oo p-2% 

e^ (l-e-^) 2 (l- 
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00 z l e~ l ^ 1 ~ m ^ i) °° z l e~ l( ~ Vl ~ V3 ^ 2 ' 1 

U *V»fe § 12l% + '" 

^(^e-fa-^/2)) pi ( ze -fa-^/2)) ^ 

77^773 12?73 

2.1.2 Oblate shape potential 

Next we discuss the case of anisotropy corresponding to an oblate shape potential U\ > uj 2 = 
U3 ( k\ = 1 < k 2 = k 3 ). In this case, u> = uj\ and 

E n = uM + u 2 (\ 2 + \ 3 ). (17) 

The number of particles excited in these dimensions are given respectively by 

No = j^-, (18) 

1 — z 

00 1z 

N t = V , (19) 



N2 = E 



"2 = 



l,n 3 =l e Z 



M=2 c ^ 



iv 3 = E 



ni=l,n 2 =0,n 3 =0 c * 



oo oo 



= V V (M + 1) " (21) 

m=i Af=o c * 

where the factor 2 in Eq. (|i"9|) is due to the symmetry between the second and third axis. 
Compared to Eq. (0), we obtain 

= 3ft(«-*) + ,.. (22) 

in the present case. 

The number of particles excited two-dimensionally on £2 — £3 plane can be written as 



n 2 = E E ( M - 1)^"' % 



Z=l Af=2 
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(23) 



For three dimensional excitations, Eq. (|2"T| ) gives 
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2.1.3 Maximally anisotropic potential 

For anisotropies cji > cj 2 > oj% with fci = 1 << k 2 « k 3 , u = uj\ and 

E n = uM + W2A2 + CU3A3. (25) 

The number of excited modes in the corresponding dimensions can be defined by 

Ni = £ " > (26) 



X' 



^2 = E r^ , (27) 

n 2 =l,n 3 =0 c * 



00 z 

n 3 = y 

ni=l,7i2=ri3=0 

^ ^ ^ O pAf?yi+A2»;2+A3»73 _ 7 ' ^ ' 

A 2 =0A 3 =0M=1 ^ c ^ 

For the two dimensional case, following Eq. (|l5l), 
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where k = k^jki- 

For three dimensional excitations, Eq. (|28|) gives 
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3 Finite Size Effects and Dimensional Crossover Be- 
havior 

3.1 Bulk behavior 

The bulk three-dimensional condensation temperature is defined in the thermodynamic limit 
?7i — > (i — 1,2,3) and N — > oo while rji/r] 3 and 772/773 fixed. The dominant term is given 
by the first term in N 3 and the critical temperature satisfies 

N= TsD C(3). (31) 

UJiUJ 2 UJ 3 

Therefore 

T3D -{-wr ■ (32) 



Two-dimensional limit is given by 772,773 "~ * 0) N — > 00, but 771 >> 1. In such a case, the 
dominant term in particle number is N 2 and the critical temperature T 2 d is defined by 

N = lA g ^ e -(^+^)/2T 2D) = ^i C ( 2 ) + . . .. (33) 

uj 2 uj 3 LO2OJ3 

Thus we have 

— p^f. (34) 

For one- dimensional limit 7/3 — > 0, A" — > 00, but 77x5772 >> 1, the dominant term in particle 
number is N\. The condensation temperature is defined by 

N = IlE gi ( e -^T 1D y (35) 

u 3 



To leading order in i] 3 1 , Eq. (p!4|) can be approximated as iVi ~ z^l°S^7- Thus the 
one-dimensional condensation temperature Ti£> is defined by 



OJ 3 <-> W3 



which gives 



AT = log , (36) 

^3 ^3 



TlD ~ M^o (37) 



for large N, Note that in the thermodynamic limit N —>■ 00 while Nu 3 fixed, T\d vanishes. 

In Figure 1, we plot Tm, T 2 o, and T 3 n as a function of the anisotropy parameter k 3 . 
In an isotropic and weakly isotropic case {k 3 ~ 1), T 3 e> < Tin,T 2 n and the condensation is 
directly into the ground state. As anisotropy is increased (k 3 » 10 3 ), T\d,T 2 d < T 3 d is 
achieved. This is the regime where various multistep behaviors can take place. 

In terms of bulk condensation temperatures we obtained in Section 3.1 as T\r> = Nuo 3 / log(2AQ, 
T 2D = (Nu 2 u 3 /2((2)) l /\ and T 3D = (A^c^/C^)) 1 / 3 , the conditions (A) T 1D « T 2D , 
(B) T 2D « T 3D , and (C) T 1D << T 3D give constraints for k 3 and k as 



(A) k » 

(B) k 3 2 /n » 



NC(2) 
(log(2iV))2 : 

A^(3) 2 

C(2) 



3 



{C) hK ^(idw- (38) 

Since (B) T 2D « T 3D implies T 3D « u>i, three-step BEC never occurs in harmonic traps. 
In Figure 2, different condensation behaviors corresponding to various anisotropy parameters 

9 



oui/u>2 and u>2/uj 3 are shown. The vertical axis corresponds to the prolate-shape potential 
studied in Section 3.2.1 (this case was studied in [0]). In such a potential, two-step BEC can 
be seen. The horizontal axis corresponds to the oblate-shape potential discussed in Section 
3.2.2, where we show that there is no multistep condensation in this case. The more general 
class of anisotropic case will be discussed in Section 3.2.3. The combined effect of dynamical 
dimensional reduction and two-step BEC in such a potential can appear in three steps. 

3.2 Dimensional crossover and condensation 

For a highly anisotropic trap, the three dimensional crossover temperature T 3D should be 
reached when the correlation length is in the order of the size of the ground state wave 
function in the most confining direction. Spreading of the wave function can be characterized 
by Li = Jh/muJi (for i = 1,2,3) |29| , p0| . Hence the above condition is equivalent to 

£(^3d) ~ ^edB/ ' VH ~ Li, where XodB = h/y2wmkT is the thermal de Broglie wavelength, 
and £3 = \T£ D — T 3D \/T 3D . This will give us the crude estimate of T^ D as 

l^-T 3D |~(^y /3 T 3D . (39) 

3.2.1 Two-step condensation 

For a prolate shape potential discussed in Section 2.1.1, we expand the whole particle spec- 
trum with respect to 771 and 773 and obtain 

N = N + Nx + N 2 + N 3 

gi (ze-^ 2 ) , 2<7 2 (ze-W2) g 3 { ze ~^) 



m mm mm 

m viva m m 

9o(z)r] 3 2g 1 (z)rj 1 g 2 (z) 
773 2 771773 2 77^773 

9o(z) fViV + 9i(z) Vi _ 9o(z) Vi 3 
77x773 V 2 / 77^773 2 77^773 6 



+ ^Mv) +?P^-^^- + ---. (40) 



This expression can be simplified to give 



9s(z) g 2 (z) gi(z ) 
V\ Vi 2 



A^3 = ^ + ^ + ^F + ---. (41) 



Writing z = e * and expanding Eq. (|41f) with respect to give 



NV3= m + m.m + ..., (42) 

Vi Vi Vi 
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where we used an asymptotic expansion of the Bose-Einstein function g${e a ) ~ £(3) — 

C(2)a + 1/2(3/2 - loga)a 2 + • • • and g 2 (e~ a ) ~ C(2) + (log a - l)a H for small a [§]. 

Correlation length £ of an ideal Bose gas is given by £ = XedB/^V^P ||31| . In terms of 
scaling parameters Xj(T) = (f>(T)/r)i(i = 1,2,3), the above argument implies that Tg*^ is 
achieved when xi(T^ D ) = c\, where C\ is some constant in the order of unity. Inserting this 
into Eq. (f42|) gives 

jy = Jg). + C(2)(l-ci) + ... at T = I ~ (43) 

Thus we have 

B» = 1 + *-i eg) fM 1/3 (44) 

r 3D + 3 c(3) 2 ' 3 W 

This result gives the same correction term proportional to [k 3 /N) as in Eq. fl39|) obtained 
by heuristic arguments. 

The one-dimensional condensation temperature Tn> is defined in Eq. ([35]) in the limit of 
small 773 and the vanishing reduced chemical potential 0; finite size effects on Tm originate 
in finiteness of both 7/3 and <fi. At one-dimensional crossover temperature T* D , correlation 
length reaches the size of the ground state wave function in the least confining direction, 
namely, along the third axis in the present case. Thus (,(T* D ) ~ L 3 or equivalently x 3 (T* D ) = 
C3 = 0(1). Then from the second line in Eq. (pEDp, we obtain 

N = ^(e-d+^^/a) | 2g 2 {e-«* + ***»l 2 ) | g 3 (e-^ + ^) | ^ 

V3 77i773 T/^7/3 

In the limit 7/3 —» 0, only the first term dominates and we obtain the crossover temperature 

as 

UJ 3 

^log . ±Lm_ . (46) 



For large N, this gives 



u 3 (1 + 2c 3 V 3 



1D Iog[2JV/(l + 2c 3 )] ■ l l) 



Note that T? D > T 1D holds. 

The conditions (C) in Eq. (RBI) and uj\ < Tm are satisfied if 



A < ^ < TT^m- (48) 



(log(2iV)) 3 / 2 d (log(2iV)) 
11 



In such a case, two-step condensation leading to the condensation into the ground state 
can be seen whereas the system is effectively still three-dimensional. In Figure 3-6, the 
condensation fractions Ni/N{i = 0, 1, 2, 3) as a function of temperature are plotted. At high 
temperature, three-dimensionally excited states dominate, as expected from the density of 
states which grows as M 2 , where M is the number of degeneracy of an isotropic harmonic 
oscillator appeared in Eq. @. 

In the isotropic case (Figure 3), condensation is only into the ground state. Due to the 
finite size effects, condensation already starts before the critical temperature is reached. In a 
strongly anisotropic case, as in Figure 4, two-step condensation can be seen. T 3 * D determines 
the onset of condensation into one-dimensionally excited states. At T 3 * D the ground state 
fraction is negligiblly small. Condensation into the ground state will not start until T* D is 
reached. 

In the multistep process peculiar to the highly anisotropic system, when the correla- 
tion length reaches the size of the system, the dynamics shows the crossover to the low- 
dimensional one before the actual phase transition occurs. In this sense, the critical temper- 
ature is never observed in such a process and the directly relevant quantity to the observation 
is the crossover temperature, the temperature at which the finite size correction sets in. For 
practical purposes, this is often replaced by including the finite size correction as the term 
proportional to the power of 1/iV whereas the chemical potential is set to the ground state 
energy [fffil . Strictly speaking, however, since the chemical potential never reaches the ground 
state energy in the finite system, the meaning of this correction has some ambiguity. The dif- 
ference between the crossover temperature and the finite size-corrected critical temperature 
in the present axially symmetric trap case is given by 

AT = ^(2WM 1/3 m 

T 3D ' 3 C(3) 2 / 3 \Nj • { } 

While this is fairly small for an isotropic or weakly anisotropic case: AT/T 3 d ~ 0.024 for 
ks = 1, it is no longer so for a strongly anisotropic case: AT/T^o ~ 0.24 for k% = 10 3 as 
used in Figure 4, where c\ = 1, N = 10 4 for both cases. The ordinary finite size correction 
significantly underestimates the results in the latter case. 2 For these reasons, we focus our 
discussions on crossover temperatures in the present work. We should also note that there 
is a slight amount of ambiguity in the choice of Ci(i = 1,2,3). In general, the correlation 
length £ is a complicated function of the temperature away from the critical value and the 
reliable choice is obtained by the numerical fitting. We will simply put q(i = 1, 2, 3) = 1 for 
our comparison with numerical data for brevity. 

3.2.2 Two-dimensional condensation 

For an oblate shape potential discussed in Section 2.1.2, assuming rjx » 1 and we obtain 

N = N + N x + N 2 + N 3 

2 Since the boundary condition in the harmonic potential corresponds to the one in the Neumann boundary 
condition, the surface correction increases the density of states and hence decreases the critical temperature. 

12 
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Two-dimensional crossover temperature T 2D can be defined when the correlation length 
reaches the size of the ground state wave function along the second axis. Thus we have 
£(T 2D ) ~ L 2 or equivalently x 2 {T 2D ) — c 2 = 0(1). Then from Eq. (|50|) we obtain 
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92{e 
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at T = T 2 * D . Expanding in terms of 772 and we obtain T 2 * D as 

C(2) 1 



(51) 
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V2 



V2 
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for 772, r/3 < < 1 . For large N and c 2 = 1 , this gives 



+ 21og(^) 
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+ 



(52) 



it 



1/2 



2D 



r : 



1 + 



2D 



M( 2 ), 



log 



iV 



X(2), 



(53) 



As explained in Section 3.1, condensation into N 2 does not occur in harmonic traps. 
For an oblate shape potential, the system dynamics freezes out along the first axis at T = 
u)\. Therefore the dynamics of the system at T < uj\ is two-dimensional. Ordinary two- 
dimensional BEC can still be observed as long as T 2 £> < u>±. This condition requires k 3 > 
(iV/C(2)) 1//2 . In Figure 5, two-dimensional BEC in this parameter regime is shown. 



3.2.3 Three-step dimensional reduction 

For a general class of anisotropic potential discussed in Section 2.1.3, three-step process can 
be observed. We expand each iV, with respect to r]i,r] 2 ,r] 3 and obtain 



-m/2\ 



go{z) + 



gi(ze 



+ 



92{ze 



-'12 



/2 ) 



K 



V3 W3 2 I 

g x (ze-^-*n-m)/2} 2 + %2 + %2 



9o{ze 



-mi 



/2-, + 93[ze_ 



-(m-m-va)/ 2 ^ 



^1^2% 



24 



vimm 



+ 



&»(*J M^J 92\z) rn-r] 2 -ri 3 



riimm ^ 2 r/ 3 r/ir/2% 



(54) 
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Expanding Eq. (|54| ) with respect to 

iV = 



C(3) | C(g) 



- log z gives 

1 1 1 

+ + 



m^m 



V1V2 mriz 773771 



t?i 772^3 



(55) 



As discussed in Section 3.2.1, xi(T 3 * D ) = c\ = 0(1) holds at T = T^ D . Inserting this into Eq. 
(|55|) gives 



N 



C(3) + C(2) 



77i772?73 



1 I /l 

+ 



771772 772773 V2 



Cl 



1 



773771 



+ ••• at T = T* . 



(56) 



Thus we have 



T JD _ 1 , Cl - 1/2 C(2) ffeAfc 



1/3 



7? 



3D 



C(3) 2 / 3 V N 



(57) 



The crossover temperature T^ D associated with two-dimensional BEC can be defined if 
T 2 * D < cjx.0 Then from Eq. (|54D we obtain 



iV 



M e 



-772(1/2+02)^ 



+ 



*- + 

7/2773 

log(l - e-^l 2 )^ 
247/27/3 



#i(e" 



-r? 2 (l/2+C2)^ 



-m/2 



% 



+ 



771772773 



+ 



(5? 



at T = T^zj- When rji becomes large, only the first two terms dominate and we obtain the 
crossover temperature as 



N 



C(2) 1/2 + c 2 | (l/2 + C2)log[r7 2 (l/2 + c 2 )] log(77 2 c 2 ) | 



772773 773 

T? D ((2) T1 D 
1^2^ 3 ^3 

+ (1/2 + c 2 ) log 



V3 



>/3 



[(1/2 + Pa) 



T 



27J 



+ log 



^2(1/2 + 02)/ & Vo;2C2 
for 772 , 773 < < 1 . For large N and c% = 1 , we simply have 



+ • 



(59) 



A 2D 



H 



1/2 



T; 



27J 



^C(2), 



log 



iV 



X(2)/ ' 



(60) 



One-dimensional crossover temperature T* D has the same form as in the prolate shape po- 
tential case. 



Thus the system behaves effectively two-dimensional (EIRD=2) at T = T^ D . 
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The conditions (A) and (B) in Eq. ( [38]) give the constraints for anisotropy parameters 
for three-step behavior to be observed: 



K > 



(log(2A)) 



h > iw' (61) 

If T 2 d « Tzd holds, since it also implies T 2 d « W\, excitations along the first axis will 
be dynamically suppressed at T < u>i making the system effectively two-dimensional before 
BEC sets in. Furthermore, if T 1D « T 2D is also satisfied, BEC occurs in two-steps, one at 
T 2 d and the other at Tn>. In Figure 6, above senario of three-step dimensional reduction 
is numerically realized. Three-dimensionally excited modes dominant in higher temperature 
are dynamically suppressed at T < uj\ followed by two-step BEC. When the condensation 
into the ground state sets in, the effective dimension of the system is still two. From the 
form of critical temperatures in Eqs. (|44]) and (|53|), we see that high anisotropy and the 
small number of atoms have similar consequences. 

It is useful to see how the weak interaction effect modifies the above arguments.^] The 
shift of the critical temperature due to the interaction effect is given in |33[ as AT^/T 3 d ~ 



aN 1 / 6 /L, where a is the s-wave scattering length, L = Jfi/mVt. The shift due to finite 

size is given in Eq.(|5"T|) as AT 3 * D /T 3 £) ~ (k 3 /N) ' . For a reasonable choice of parameters, 
N = 5000, k 3 = 10, a/L = 0.001, we obtain AT^/T 3D ~ 0.003 and AT* D /T 3D ~ 0.13. 
Thus the interaction effect is negligible compared to the finite size correction. 

Comparing the zero point energy in the most confining direction with the interaction 
energy, if u)\ > n U, where n is the ground state density and U = <iiifi 2 a/m, the criteria for 
the system to behave as effectively two-dimensional is still given by u)\ > T. For the choice 
of value n U/kB ~ 110 [nK] for a sodium atom, this gives uj\ > n^U ~ 10 4 [1/sec]. We 
should note that the collision effect even in such a case cannot be strictly two-dimensional 



34jl ; for high energy atoms, the interaction vertex has the form of the three-dimensional 
one. Nevertheless, these effects will be still suppressed for a « L\ and a small num- 
ber of particles whereas the finite size correction behaves as AT 2D /T 2 £> ~ {n/N) and 
AT^ d /Tid ~ 1/logiV. Therefore the effect of interaction to the particle number and the 
critical temperature remains small and the system is still expected to show the multistep be- 
havior under these conditions. At T < u>i, the most notable difference would be the absence 
of long range order due to interactions in this effectively two-dimensional system. In this 
case, the system can show quasi-condensation, the condensation with the fluctuating phase 
at near T 2 £>. The true condensation with the constant phase will be acheived at lower tem- 
perature. One example of quasi two-dimensional system, a gas of spin-polarized hydrogen in 



liquid helium, is known to exhibit Kosterlitz-Thouless transition [35, 36]. Manifestation of 



the crossover from BEC to Kosterlitz-Thouless transition [37| during the multistep process is 



of particular interest to study. But even in the presence of such a transition, we expect that 



We recover ordinary units in this discussion for quantitative comparison with other literatures. 
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the particle occupation number which is insensitive to the phase information still behaves 
similar to the ideal gas. The interaction contributions used above are calculated within the 
framework of the mean field theory. Thus we conclude, apart from the critical regime where 
the mean field theory fails, interactions do not essentially alter the multistep process for a 
small number of atoms in highly anisotropic traps. 

From the arguments given in this section, the sum of most relevant terms in Eq. fl54|) 
around each crossover temperature will give the simple expression of the total number of 
atoms as 

N = N + JVi + iV 2 + N 3 

„ go(e -») + ME!) + *<£^ + i^l . (62) 

V3 Ws r niV2V3 

Making use of the fact that <fi varies as a nonvanishing function of the temperature for 
the finite system, we define N 3 (X) = g 3 {e x ^)/f]ii] 2 ri 3 . The relation for the Bose function 
^n(e~*) = -g n -i(e~' t> ) @j allows us to write 



N 2 {\) = 

iVi(A) = 



g 2 (e^) _ 1 dN 3 (X) 

r] 2 r] 3 xi dX 
gi (e*) _ 1 d 2 N 3 (X) 



r] 3 xix 2 d\ 



2 



Thus the total number of atoms is given by 

N = ^^!^ + ^^zl) + i^) +JV3( _ 1) . (64) 

X\X 2 x 3 dX 6 X\X 2 dX z x\ dX 

At T ~ T 3 * D , xi ~ 1 and a?2, ^3 >> 1, then we have 

iV ~ Ajv 3 (-1) + jV 3 (-1). (65) 

At T ~ T 2 * D , x 2 ~ 1, xi << 1, and x 3 » 1, we have 

N „ — JV 2 (-1) + JV 2 (-1). (66) 

uA 

At T ~ T* D , X3 ~ 1 and xi, x 2 << 1, we have 

AT ~ ^.jViC-lJ + iViC-l). (67) 

aA 

These results suggest that each condensation fraction Ni behaves similarly as a function of 
the temperature in the neiborhood of its characteristic temperature T* D . We simply write 
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this fact as Ni(T) ~ F(T/T* D ) for i = 1, 2, 3, where F is a function independent of i. It also 
implies that each component has a similar shape in the logarithmic T scale. The occupation 
number of each dimension is plotted in the logarithmic T scale in Figure 7. Note that this 
derivation relies on the special property of the Bose function which determines the density 
of states of an ideal gas trapped in the harmonic potential. Whether the same result holds 
or not in other systems with different density of states is not obvious. 

In conclusion, finite size effects on the Bose-Einstein condensation of an ideal gas in a 
strongly anisotropic trap give rise to various different types of multistep behavior depending 
on the degree of anisotropy. In an isotropic trap, BEC into the ground state always begins 
while the system is effectively three-dimensional, i.e. T^d > i>j\ = uo 2 = 003. In an anisotropic 
trap, in addition to the BEC which may occur in multisteps, EIRD will also decrease in steps 
as the temperature is lowered. The combined effect of these leads to the appearant multistep 
behavior. The existence of the intermediate condensation into one- dimensional space can be 
traced back to the logarithmic divergence of the one-dimensional occupation number in Eq. 
(J36|). This means that, when the trap is loosened in one direction, the particles tend to occupy 
quantum states along this direction with more likelihood than along other directions. Thus 
one-dimensionally excited modes in this direction will dominate multidimensional excitations 
spread in other directions (Aq » N 2 , A/3) and the thermodynamic behavior of such a system 
is characterized by T\d even though effective dimension of the system is still three. Note 
that the same mechanism is responsible for the nonexistence of BEC in one-dimensional 
harmonic trap in the ordinary thermodynamic limit. For the same reason, the intermediate 
condensation into two-dimensionally excited modes can be observed in a rectangular cavity 
PTDj , ^| , where T 2 d does not exist in the naive thermodynamic limit. Three-step BEC can 
take place only in such a system. Away from the thermodynamic limit, the temperature 
dependence of the chemical potential around T 1£) , T 2D , and T 3D causes similar crossover 
behaviors in condensation fractions Aq, A r 2 , and A3 as a function of the reduced temperature. 

Atom trap experiments probing the two-step BEC are realizable in Ioffe-Pritchard type 



magnetic traps or in optical dipole traps ||24j| . The similar type of device can be used to study 
multistep behavior discussed here, although it is difficult to achieve BEC in a maximally 
anisotropic trap with our current cooling technique. Further progress in a trapping device 
may be required. The basic mechanism of multistep dimensional crossovers discussed here 
can be applied to many other bosonic systems and should be amenable to future experimental 
verification. Quasi low-dimensional systems realized in the optical lattice or waveguide are 



the promising option for testing such processes [17|, |18fl . Also of interest is the kinetics of 
multistep behavior where the correlation length and the thermal de Broglie wave length are 
related in a nontrivial manner. This is currently under investigation. 
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A Dynamical Symmetry in Anisotropic Harmonic Os- 
cillator 

It is known that dynamical symmetry of anisotropic harmonic oscillators has the reducible 
representation which can be written as the multiple of irreducible representation of SU(3) 
symmetry, the symmetry of isotropic harmonic oscillators. 

as 



For a given A = (Ai, A2, A3), we can define a set of many boson annihilators J38 



4 A = ^y[|l^T^ fOT « = 1,2,3 , (Al) 

where hi = ajai are boson number operators. Many boson annihilators and creators satisfy 
the following commutation relations 






S^. (A2) 



Rewriting the Hamiltonian M) as 



H = ^E(A A , 4 f } -\jl<h - 2\ - 1), (A3) 

the corresponding energy eigenvalue becomes 

E n = lu(M+ ~) - -£>(** - 2A * - !)■ ( A4 ) 

2 2 »=i 

This gives the alternative method to derive energy eigenvalues given in Eq. (^). 

Thus the reducible representation (^) of the original Hamiltonian leads to the cluster 
of isotropic harmonic oscillators. With this decomposition, it is possible to understand the 
mechanism of condensation in anisotropic traps in terms of condensation in isotropic traps 



B Mathematical Formulas 

By taking derivatives of 

00 
E e~ MTI = T^—r, ( B1 ) 

by i] on both sides, we obtain 

y Me- Mr > = - C — - (B2) 

5-L, 1-e-" 2 v ; 



1 — p -7 ? 



e"" 



M=l 



and 



Af=l V- 1 e / 



For a small 77, we have an expansion 



1 - e-* 2 sinh(r//2) 77 24 

Similarly, we have 



■n/2 p -v/2 np-ri/2 

VC +•••. (B4) 



-k V -{k-l) V /2 1 -{k-l)r,/2 

-( k + r)^r, — + ■•• ( B5 ) 



and 



(1 - e- kr >)(l - e- T >) krf v k' 24 



-m P -(vi-V2-V3)/2 



[I - e-»B)(l - e-%) 8sinh(r7!/2) sinh(?] 2 /2) sinh(77 3 /2) 

e -(m-V2-m)/ 2 e -(m-m-V3)/2 ^2_j_2_|_" 



^7i^2?73 24 ?7i^2^3 



+ ■ • • . (B6) 



From Eq. 



00 -Ip-lr] 00 Z p-lri/2 1 

E^ - g £ V = ~X«(l-«^)- (B7) 
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Figure Captions 

Figure 1 The condensation temperatures T 1D (solid curve), T 2D (dashed curve), and T 3D 
(dot-dashed curve) are shown as a function of k%. N = 10 4 , u)\ = 0.5, k 2 is also varied as 
k 2 = k 3 . The logarithmic scale is used for both axis. 

Figure 2 Different condensation behaviors corresponding to different anisotropy parameters 
uj\juj 2 and io 2 jiOj > are indicated. N = 10 4 is chosen. The logarithmic scale is used for both 
horizontal and vertical axis. ((2) and £(3) are approximated as one for simplicity. 

Figure 3 The condensation fractions Nq/N (solid curve), N\/N (dashed curve), N 2 /N (dot- 
dashed curve), N 3 /N (dotted curve) as a function of the temperature for an isotropic trap 
are plotted. The same symbols are used in Figures 4-7. uj\ = 0.1, u 2 = 0.1, U3 = 0.1, and 
iV = 1000 are chosen. T c = 0.94 is the three-dimensional critical temperature in Eq. 



Figure 4 The condensation fractions for a prolate shape anisotropic trap are plotted. u>i = 
0.3, uj 2 = 0.3, u>3 = 0.0003, and N = 10 4 are chosen. Two-step condensation occurs in this 
case. T c = T 3 * D = 0.61 from Eq. ( |44"|) and T* D = 0.34 from Eq. (|47D . c\ and C3 are set equal 
to one. Note that EIRD=3 during the whole two-step BEC process. 

Figure 5 The condensation fractions for an oblate shape anisotropic trap are plotted. U\ = 
0.3, u 2 = 0.002, lu 3 = 0.002, and iV = 10 3 are chosen. T c = T 2 * D = 0.057 is the two- 
dimensional crossover temperature in Eq. ( J53|) . c 2 is set equal to one. The system shows the 
dimensional reduction from three to two dimension at T = uo\ = 0.3. Hence the accumulation 
of particles into N 2 from A^ 3 is not the result of condensation. At T < u>i, the system show 
the ordinary two-dimensional condensation into the ground state. 

Figure 6 The condensation fractions for a maximally anisotropic trap are plotted. Three- 
step dimensional reduction can be seen. uj\ = 0.3, u 2 = 0.02, L03 = 0.0004, and JV = 5x 10 3 
are used. T 3 * D = 0.17, T 2D = 0.06, and T* D = 0.02 are crossover temperatures defined in 
Eqs. (|57|) , (|60|) , and (|47|) , respectively. c\ = c 2 = C3 = 1 are used. The system shows the 
dimensional reduction at T = u\ = 0.3 similar to Figure 5, and behaves two-dimensionally 
near T < 0.3. Therefore, the accumulation of particles into N 2 from iV^ is not the result of 
condensation. For T < u>i, the system shows two-step BEC at T 2D into one-dimensionally 
excited states and at T* D into the ground state. 

Figure 7 The logarithmic T scale is used for the three-step behavior shown in Figure 6. 
The same parameters as in Figure 6 are used. 
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